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The dispersion properties of rectangular metallic waveguides periodically loaded by uniaxial res- 
onant scatterers are studied with help of an analytical theory based on the local field approach, the 
dipole approximation and the method of images. The cases of both magnetic and electric uniaxial 
scatterers with both longitudinal and transverse orientations with respect to the waveguide axis are 
considered. It is shown that in all considered cases waveguides support propagating modes below 
cutoff of the hollow waveguide within some frequency bands near the resonant frequency of the 
individual scatterers. The modes are forward ones except the case of transversely oriented magnetic 
scatterers when the mode turns out to be backward. The described effects can be applied for the 
miniaturization of the guiding structures. 
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I. INTRODUCTION 

Recently, a very unusual waveguide was proposed by R. 
Marques in [1] and then extensively studied by S. Hrabar 
in [2]. It is a rectangular metallic waveguide periodi- 
cally loaded by resonant magnetic scatterers, so-called 
split-ring-resonators (SRR:s) [3, 4], which are also used 
as components of a realization of the left-handed medium 
(LHM) [5], composite with negative permittivity and per- 
meability [6, 7]. The geometry of the Marques waveguide 
(MW) is presented in Fig.l. The SRR:s in MW are ori- 
ented so that their magnetic moments are orthogonal to 
the waveguide axis and to one of the walls. 




FIG. 1: Geometry of a subwavelength Split-Ring-Resonator- 
loaded metallic waveguide 

The MW support a propagating mode within a fre- 
quency band near the resonance of SRR:s even if it is 
located below cutoff frequency of the hollow waveguide 
[1, 2]. The transversal dimensions of the waveguide hap- 
pen to be much smaller than the wavelength in free space. 
Thus, loading by SRR:s makes waveguide subwavelength 
and provide unique method for miniaturization of guid- 
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ing structures. The mode of MW is backward wave (the 
group velocity is negative) . This effect was interpreted in 
[1] in terms of the effective LHM to which such the loaded 
waveguide is apparently equivalent. The empty waveg- 
uide was considered as an artificial electric plasma with 
negative permittivity and the array of magnetic scatter- 
ers as a magnetic material with negative permeability. 
This interpretation is not completely adequate because 
it cannot explain why the effect disappears in the case of 
loading by isotropic magnetic with negative permeabil- 
ity. Really, it is clear that the hollow waveguide filled by 
isotropic magnetic with negative permeability does not 
support guiding modes. Meanwhile, the doubly-negative 
medium having the isotropic negative permittivity and 
permeability would support propagating backward waves 
[8]. Also, the LHM interpretation is not instructive in our 
opinion since it does not allow to notice possibilities to 
obtain propagation below the cutoff frequency of hollow 
waveguide with help of the other loadings that SRR:s. 

The goal of the present study is to give an adequate 
explanation for the extraordinary propagation effect in 
MW and to suggest another loadings which would lead 
to the similar effects. In this paper it is shown that the 
propagation below the cutoff frequency of hollow waveg- 
uide can be achieved with magnetic scatterers oriented 
longitudinally with respect to waveguide axis, as well as 
with electric scatterers oriented either longitudinally or 
transversally. This demonstrates that the miniaturiza- 
tion of the rectangular waveguide at a fixed frequency us- 
ing loading by the resonant scatterers is not restricted by 
the case when the scatterers are magnetic and transver- 
sally oriented. The miniaturization is possible with help 
of either magnetic or electric resonant scatterers with ei- 
ther transversal or longitudinal orientation with respect 
to the waveguide axis. Of course, the miniaturization 
can be also reached using periodically located capacitive 
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posts, however the loading by resonant scatterers is a 
qualitatively different effect. In [2] it was pointed out 
that the miniaturization obtained in this way refers also 
to the longitudinal size of the waveguide since the period 
of the loads in incomparably smaller than the wavelength 
in free space, unlike the propagation in a capacitively 
loaded waveguide, where the period of the posts is of the 
order of A/ 2. 

The mini pass band below the cutoff frequency of 
a rectangular waveguide loaded by resonant scatter- 
ers is caused by the properties of the periodical one- 
dimensional array (chain) of resonant dipoles and has 
nothing to do with doubly-negative media. The back- 
ward wave appears in a special case of transverse orien- 
tation of magnetic scatterers and is not a necessary at- 
tribute of such mini-band. It is known that a chain of the 
resonant scatterers in a homogeneous matrix supports 
guided modes. In the optical frequency range it refers to 
chains of metallic nanoparticlcs [9-12]. At microwaves it 
refers to so-called magneto-inductive waveguides (chains 
of SRR:s) [13-16] or to chains of inductively loaded elec- 
tric dipoles [17]. The metallic walls of loaded waveguide 
perturb the dispersion of the guided mode in a chain 
but do not cancel the propagation. This is because the 
wavelength of the guided mode in the chain of resonant 
scatterers is dramatically shortened compared to that in 
the matrix. As a result, the energy of a guided mode 
is concentrated in a narrow domain around the chain, 
and the interaction between the chain and the waveguide 
walls turns out to be not critical for the existence of the 
guided mode. 

The paper is organized as follows. In the Section II the 
dispersion properties of the chains of resonant scatterers 
located in free space are considered. The known results, 
obtained in the work [12] by numerical simulations, are 
reproduced with help of analytical theory based on lo- 
cal field method. This is a necessary part of the work 
in the view of comparison with the case of the loaded 
waveguide. The coincidence with known results can be 
considered as a validation of our approach. In the Section 
III the dispersion properties of the chains located inside 
of the rectangular waveguide are considered using two 
approaches: an accurate method of local field (as in Sec- 
tion II) and an effective medium filling approximation. 
The Section IV is devoted to comparison between prop- 
erties of the chains and loaded waveguides. The Section 
V contains concluding remarks. The details of the local 
field theory are given in Appendix. 

In the present paper we consider both magnetic and 
electric resonant uniaxial scatterers. As an example of a 
magnetic scatterer we have chosen the SRR:s [3, 5, 18] 
(see Fig. 2. a). The electric dipoles are represented in our 
work by the short inductively loaded wires (ILW) [19] (see 
Fig. 2.b). Any individual scatterer can be characterized 
by polarizability relating the dipole moment (magnetic or 
electric) with the local field (magnetic or electric external 
field acting to the scatterer) . This polarizability is scalar 
since the only possible direction of the induced dipole 



moment is possible for a scatterer with given orientation. 
The details concerning calculation of the polarizabilities 
for SRR:s and ILW:s are presented in Appendix A. 




21 



FIG. 2: Geometries of resonant scatterers: a) Split-Ring- 
Resonator, b) inductively loaded wire dipole. 

The inverse values of the polarizabilities a{u>) and 
a e (uo) (see Appendix A, formulae (Al) and (A4)) of 
SRR:s and ILW:s have the same dependencies on fre- 
quency within the resonant band: 

Re{a" V)} = A- 1 (ji - l) . (1) 

Here A is amplitude and uj is resonant frequency, the 
parameters determined by the geometry of the scat- 
terer. Notice, that the result (1) is also valid for a silver 
nanosphere in the vicinity of its plasmon resonance [12]. 
Thus, it is clear that there is no principal difference be- 
tween the dispersion properties of the chain of SRR:s and 
ILW:s (at microwaves) or silver nanospheres (in the op- 
tical range). 

II. CHAINS OF UNIAXIAL RESONANT 
SCATTERERS 

Let us study dispersion properties of the linear chains 
with period a formed by resonant scatterers. We will con- 
sider only two typical orientations of scatterers: longitu- 
dinal and transverse. The geometries of the structures 
are presented in Fig. 3. The case of longitudinal ori- 
entation was analyzed in [17], and the both longitudinal 
and transverse orientations were studied in [12]. In the 
present section we reproduce the main results of these 
works with help of the local field approach. 

A. Basic theory 

Without loss of generality we can restrict consideration 
by the case of the chain of magnetic scatterers (SRR:s). 
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FIG. 3: Chains of resonant scatterers. a) Longitudinal orien- 
tation, b) Transverse orientation. 



It allows to get dispersion equation for the chains under 
consideration: 

/ioa(w)^ 1 = C d (uj, q, a), (7) 

where 

C d (u,q,a) = G dd (amx Q )e-^ am . 

In the Appendix B we provide expressions (Bl) and (B2) 
which we use for effective numerical calculations of inter- 
action constants C x and C y , corresponding to transverse 
and longitudinal orientations of scatterers, respectively. 



The chains of electric scatterers are dual structures to the 
considered ones and have completely the same dispersion 
properties. 

The spatial distribution of dipole moments of SRR:s 
corresponding to an eigenmode of a chain is determined 
by a propagation constant q: M n = Me" M ™. Following 
to the local field approach the dipole moment M of a 
reference (zeroth) scatterer can be expressed in terms 
of the magnetic field Hi oc . acting to it: M = aHf oc , 
where Hf oc = (Hi oc . • d) is the projection of the field on 
the direction of the scatterer (d = Xo for longitudinal 
orientation of scatterers and d = yo for transverse one). 
This local field is a sum of partial magnetic fields H,„ 
produced at the coordinate origin by all other scatterers 
with indexes m^O: Hi oc . = H m . 

The magnetic field produced by a single scatterer with 
dipole moment M m at a point with radius vector R is 

given by dyadic Green's function G(R): 

H m (R) = Mo l W(R)M m , (2) 

where 

—jkR 

^ R ) = ( fc2T+vv )^F- ^ 

Since all dipole moments of the chain are oriented along 
d it is enough to use only the dd component of dyadic 
Green's function. So, we replace (2) by the scalar expres- 
sion: 

H^(R)= At o 1 G dd (R)M m , (4) 

where 

/ fp \ e -jkR 

and d means x for the longitudinal case and y for the 
transverse case. 

Finally we obtain the expression for the field acting to 
the reference scatterer in the form: 

Kc. = E Gdd(amx )e~™ am M. (6) 

m/0 



B. Analysis of dispersion properties 

The dispersion diagram for guided modes can be ob- 
tained solving transcendental dispersion equation (7) 
with interaction constants given by expressions (Bl) 
and (B2). Geometrically, dispersion curves correspond 
to the lines of level where the surface plot of func- 
tion Re{C X}V (uj, q)} is crossed by /^oRe{a" 1 (cj)}. Note, 
that only the solutions with k < q < lixja — k for 
k < it /a correspond to guided modes: For \q\ < k, 
Ivci{C Xt y(u>, q) — Hoar } 7^ and dispersion equation (7) 
has complex solutions corresponding to leaky modes (see 
details in Appendix B). 
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FIG. 4: Dependence of Ke{C x a 3 } on normalized frequency 
ka/n and propagation factor qa/n 

The dependencies of C x and C y on on normalized fre- 
quency ka/w and propagation factor qaj-n are shown 
in Fig. 4 and 5, respectively. The interaction con- 
stants vary in [—1, 1] range except the case of C v with 
q close to k, which has logarithmic singularity at the 
light line q — k. The function /ioRc{a _1 (w)} decreases 
very rapidly within [—1, 1] range of values near the res- 
onant frequency ujq. It means, that guided modes ex- 
ist only within a narrow bands near the resonant fre- 
quency of scatterers. The behavior of dispersion curves 
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FIG. 5: Dependence of Ke{C y a 3 } on normalized frequency 
ka/ir and propagation factor qa/ir 



can be easily predicted from the plots in Fig. 4 and 
5. If at a fixed frequency the interaction constant de- 
cays when propagation factor increases then the disper- 
sion curve grows and the eigenmode is forward wave (the 
group velocity v g = ^ is positive), but if the inter- 
action constant grows then the dispersion curve decays 
and the eigenmode is backward wave (the group velocity 
v g = ^ is negative). From Fig. 4 it is clear that for 
any resonant frequencies satisfying to the evident condi- 
tion wo < 7r/(a^/£oMo) (corresponding to the propagation 
below cutoff of the hollow waveguide ) the longitudinal 
mode is forward wave because C x decays versus q. In the 
case of transverse modes the situation is different. While 
k$a < 0.57T (fco = Wo^/eoMo) a two mode regime holds. 
The interaction constant C y decays while q is close to 
k, but from a certain q it starts to grow. It means that 
one of the transverse eigenmodes is forward (with q ps k) 
and the other one is backward. If the resonant frequency 
is high enough (0.57T < k$a < ir) the two mode regime 
disappears and only the forward wave remains. 

The typical dispersion diagram for both longitudinal 
and transverse modes is presented in Fig. 6 for the case 
of scatterers with A — 0.1/zoa 3 and oj^a = l/^/eoMo- The 
similar dispersion diagram was obtained in [12] (see Fig. 
3 of this work) by a numerical simulation. Though in 
[12] the electric scatterers in the optical range were con- 
sidered, but we consider the magnetic scatterers in the 
microwave range, the use of duality principle, and nor- 
malized frequency ka/ir and wave vector qa/ir eliminates 
this difference. The polarizability of silver nanospheres, 
for which Fig. 3 from [12] was obtained, obeys to expres- 
sion (8) of [12] which is identical to our formula (1)). 

As it was predicted, the longitudinal mode is forward 
wave and there is a two mode regime for transverse 
modes. The dispersion curve for transverse waves has 
the asymptote q = k and both leaky (q < k) and guided 
(q > k) modes exist at very low frequencies where they 
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FIG. 6: Dispersion diagram for chains of resonant scatterers: 
transverse (thick line) and longitudinal (thin line) orientations 



have almost equal wave vectors. This fact indicates the 
good matching between the radiated wave and the guided 
mode. Within the band 0.995 < ka < 1 there are two 
guiding modes at every frequency. The solution corre- 
sponding to the backward wave is close to the Bragg's 
mode (qa « 7r) whose group velocity is close to zero. The 
field of this mode is concentrated near the chain within 
the spatial region r = \J z 2 4- y 2 < 1/ \Jq 2 — k 2 ~ a. 
The same concerns the longitudinal mode within the 
band 1.015 < ka < 1.020. If the period of the chain 
is much smaller than wavelength in free space the waveg- 
uide is sub-wavelength (the field of the guided mode is 
concentrated within a cylindrical domain whose diame- 
ter is much smaller than A). Thus, the chains of resonant 
scatterers (electric or magnetic, it does not matter) form 
sub- wavelength waveguides which can support either for- 
ward or backward waves [9-17]. 



III. LOADED WAVEGUIDES 

A. Basic theory 

Let us study the eigenmodes of the rectangular metal- 
lic waveguides periodically loaded by resonant uniaxial 
scatterers. Such structures can be effectively considered 
as linear chains located inside of the metallic waveguides. 
The geometries of the four waveguides considered in the 
present paper are shown in Fig. 1 (left sides of subplots). 
The chains with period c along the waveguide axis are 
located at the center of rectangular metallic waveguides 
with dimensions a x b. The structures in Figs, l.a-d differ 
by orientation of scatterers (longitudinal or transverse) 
and their type (electric or magnetic). The first structure 
(with transversely oriented magnetic scatterers) is the 
sub- wavelength waveguide (see Fig. 1) suggested by R. 
Marques [1, 2]. The other ones are considered in order to 
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show three other possible ways to obtain miniaturization 
of rectangular waveguides. 

Note, that the chains of electric scatterers are no more 
dual to the chains of magnetic scatterers (in contrast to 
the chains in free space) due to the different interaction 
of electric and magnetic dipolcs with metallic walls. 
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FIG. 7: Transformation of the waveguide problem to the lat- 
tice one with the use of the image principle 

The dispersion equation for the chains keeps the same 

form as (7), but free-space dyadic Green's function G(R) 
(3) should be replaced by Green's function of the waveg- 
uide which takes into account the metallic walls. This 



Green's function can be determined with help of the 
image principle. This approach transforms the eigen- 
mode problem for the loaded waveguide to the problem 
of the eigenwave propagation in a three-dimensional elec- 
tromagnetic lattice formed by same scatterers. The de- 
tails of the transformation are illustrated by Fig. 7 (right 
parts of all subplots). The electromagnetic crystals ob- 
tained in such a way have orthorhombic elementary cell 
axbxc and their dispersion properties was studied in [20] 
using local field approach. Thus, we can apply the the- 
ory of the electromagnetic interaction in dipole crystals 
presented in [20] in order to study dispersion properties 
of waveguides under consideration. 

In the coordinate system associated to the axes of the 
crystal the center of a scatterer with indexes (m, n, I) 
has coordinates 'R. m . n ,i = (am, bn, cl) T . From Fig. 7 it is 
clear, that distribution of dipole moments in the lattice 
has the following form: 



{-l) m Me- jqcl 



x , 



for the case of transverse orientation of magnetic scatter- 
ers (Fig. 7. a); 



Me- jqcl 



for the case of longitudinal orientation of magnetic scat- 
terers (Fig. 7.b); 

Pm,„,J = (-l) n Pe-^x , 

for the case of transverse orientation of electric scatterers 
(Fig. 7.c); and 

P m ,n,i = (-l) m+n Pe-^ cl z , 

for the case of longitudinal orientation of electric scatter- 
ers (Fig. 7.d). 

Any of these distributions can be rewritten in terms of 
a wavevector q as e~ ■i(i" R "».'».') ) where the wavevector q 
for four cases considered above has the form (tr/a, 0, q) T , 
(0,0,q) T , (0,7r/6, q) T and (ir/a,n/b,q) T , respectively. 
This notation finally makes clear that the waveguide dis- 
persion problems reduce to those of the three-dimensional 
lattices in the special cases of certain propagation direc- 
tions. 

The dispersion equation for three-dimensional electro- 
magnetic crystal formed by magnetic scatterers oriented 
along a;-axis has the form [20]: 



-V)-C(fc,q) = 0, 



(8) 



where 



C(k,q,a,b,c) = G(R ro ,„,0e- J '^ OTO+ ^ 6n+g ^. 

(m,n, 1)^(0, 0,0) 

(9) 

We call C{k, q, a, b, c) as the dynamic interaction con- 
stant of the lattice using the analogy with the classical in- 
teraction constant from the theory of artificial dielectrics 
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and magnetics [21]. The explicit expression for C for the 
general case was derived in [20] and it is given in Ap- 
pendix C by formula (CI). 

The dispersion equation for waveguides with trans- 
verse orientation of scatterers can be directly obtained 
from equation (8) by substitution q = (w/a,0, q) T and 
q = (0, 7r/6, q) T for magnetic and electric scatterers, re- 
spectively. Also, in the case of electric scatterers fol- 
lowing duality principle /ioa _1 (o;) should be replaced by 
Soa^ 1 ^). The similar operation for case of longitudi- 
nal orientation happens to be possible only after rota- 
tion of coordinate axes: z — > x' , x — > y 1 , y — > z' , since 
equation (8) requires scatterers to be directed along x- 
axis, but for longitudinal orientation they are directed 
along z-axis. After such manipulation substitution of 
q = (q, 0,0) T and q = (q, w/a, n/b) T (in the new coor- 
dinate axes (x',y',z')) into equation (8) provide desired 
dispersion equations for waveguides with transverse ori- 
entation of magnetic and electric scatterers, respectively. 

This way we obtain the following dispersion equations 
for all loaded waveguides under consideration: 



Una x (w) - C(k, (ir/a, 0, q) T , a, b, c) = 0, 



(10) 



for transverse orientation of magnetic scatterers; 

^a- 1 ^) - C(k, (q, 0, Of, c, a, 6) = 0, (11) 

for longitudinal orientation of magnetic scatterers; 

e Q^V) - C(k, (0, ir/b, q) T , a, b, c) = 0, (12) 

for transverse orientation of electric scatterers; and 

eoo^M - C(k, (q,ir/a,ir/b) T ,c,a,b) = 0, (13) 

for longitudinal orientation of electric scatterers. 

The obtained dispersion equations are real valued ones, 
because the imaginary parts of its components cancels 
out. It can be clearly seen from Sipe-Kroncndonk condi- 
tion (A2) and following expression proved in [20]: 

Im(C) = j£. (14) 



B. Effective medium filling approximation 

The chain of SRR:s with transverse orientation located 
in the waveguide has been interpreted in the literature as 
a piece of an uniaxial magnetic medium [1, 2]. We call 
this approach as effective medium filling approximation. 
It can be applied practically to every waveguide consid- 
ered in this paper except the case of longitudinal orien- 
tation of magnetic scatterers since the uniaxial magnetic 
model does not describe longitudinal modes. This ap- 
proach provide qualitatively acceptable results which are 
compared with exact ones in the next subsection. 

Let us start from the case of transversely oriented mag- 
netic scatterers (Fig. 7. a) and consider a chain of parallel 



uniaxial magnetic scatterers as a piece of infinite resonant 
uniaxial magnetic. The permeability of such a magnetic 
is a tensor (dyadic) of the form: 

JI = ^x x + Aio(yoyo + z z ). 

The permeability [i along the anisotropy axis x, can be 
calculated though the individual polarizability of a single 
scatterer using the Clausius-Mossotti formula: 



H = no 1 + 



1 - C 3 (a,b, c)a{Lo)/fi 



(15) 



where V = abc is a volume of the elementary cell of 
an infinite three-dimensional lattice and C' s (a, 6, c) is the 
known static interaction constant of the lattice [20, 21]. 
In the case of a simple cubic lattice a = b = c the interac- 
tion constant is equal to the classical value C s — 1/(3^). 

Notice, that we should skip the radiation losses con- 
tribution in expression (A3) while substituting into for- 
mula (15). This makes permeability purely real number 
as it should be for lossless material. This manipulation is 
based on the fact that the far-field radiation of the single 
scatterer is compensated by the electromagnetic interac- 
tion in a regular three-dimensional array, so that there 
were no radiation losses for the wave propagating in the 
lattice [20, 22]. 

The dispersion equation for the uniaxial magnetic 
medium has the following form (see e.g. [23]): 



Mo(«£ + Q 2 Z ) = M(fc 2 - (&)■ 



(16) 



To solve the waveguide dispersion problem is to solve 
the dispersion problem (16) for a special case q = 
(n/a, 0, q) T . The substitution of q = (n/a,0,q) T into 
(16) gives: 



q 2 = — 
Ho 



(- 



(17) 



For the frequencies below cutoff of the hollow waveg- 
uide the expression in the brackets of (17) is negative and 
for positive fi there is no propagation. However, if \i is 
negative (it happens in accordance with (Al) and (15) 
within a narrow frequency range near the resonance of 
the scatterers, just above the resonant frequency of the 
media) q becomes real. This mini pass-band can be lo- 
cated much lower than the cutoff frequency of the empty 
waveguide with help of reduction of resonant frequency 
of the scatterers. It is easy to see from (17) that the 
mode is backward wave, i.e. 4^ < 0. This follows from 
the basic inequality ^ > (Foster's theorem). 

So, for transverse magnetic scatterers the effective 
medium filling model gives (at least qualitatively) a cor- 
rect result. Namely, it predicts a mini-band within the 
resonance band of SRRs and the backward wave prop- 
agating within it. However, this model is not accurate. 
The reason of this inaccuracy is simple. In spite of the 
low frequency of operation (the waveguide dimensions 
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are small compared to the wavelength in free space) the 
effective magnetic medium should operate in the regime 
when its period is comparable with the wavelength in the 
effective medium because q x = n/a. Such regime for an 
electromagnetic crystal can not be described with help of 
homogenization and requires taking into account spatial 
resonances of the lattice [20]. 

In the case of transversely oriented electric scatterers 
(Fig. 7.c) the effective medium filling model implies that 
the wave propagates in a uniaxial dielectric with permit- 
tivity tensor: 

t = ^x x + £o(y yo + z z ). 

The permittivity e along the anisotropy axis x is given 
by the Clausius-Mossotti formula: 

£ ~ £o ri-c s ( fl ,MKM/e„J' (18) 

The dispersion equation for such uniaxial dielectric reads 
[23]: 

e (ql + q 2 z )=e(k 2 -q 2 x ). (19) 

Solution of waveguide dispersion problem corresponds to 
the case when q = (0, n/b, q) T : 




This mode propagates only at the frequencies when per- 
mittivity takes high positive values e > £o[7r/(fc&)] 2 . In 
our case of resonant dielectric it happens within a mini- 
band just below the resonance of the medium. It is clear 
from (20) that the mode is forward wave, i.e. 4j > 0. 
This follows from the Foster's theorem 4 s - > 0. 

In the case of longitudinally oriented electric scatter- 
ers (Fig. 7.d) the solution of the waveguide dispersion 
problem can be obtained from the dispersion equation 
(19) with q = (q,ir/a,ir/b) T (the axis were transformed 
in order to have x-axis along dipoles) in the next form: 




The mode is propagating at the frequencies when per- 
mittivity is positive and rather high, or negative. This 
mode is forward wave in the same manner as (20) since 
4^ > 0. 

C. Analysis of dispersion properties 

For numerical calculation of dispersion curves using 
(10), (11), (12), (13) and (CI) we have chosen square 
waveguides (a = b = c) loaded by scatterers with the 
same parameters which were used for studies of chains: 
loo = l/(ay/eojIo), and A = 0.1 fioa? f° r magnetic scat- 
terers, and A e = O.leoa 3 for electric ones. 
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FIG. 8: Dependence of the real part of normalized interaction 
constant C(fc,q)a 3 with q = (n/a,0,q) T (corresponding to 
transverse orientation of magnetic scatterers) on normalized 
frequency ka/n and propagation constant qa/n. 



The dependence of the real part of normalized interac- 
tion constant C(k, q) with q = (7r/a,0,g) T on the nor- 
malized frequency ka/ir and on the propagation constant 
qa/n is presented in Fig. 8. This interaction constant 
corresponds to dispersion equation (10) for transverse 
orientation of magnetic scatterers. The value of Re(C)a 3 
varies within [—2,0.5] interval while the normalized fre- 
quency ka/ir is bounded by unity (which corresponds to 
the cutoff of hollow waveguide). If a value of the nor- 
malized frequency is fixed then the real part of the in- 
teraction constant is a monotonously growing function 
of qa/ir. The dependence of the real part of the inter- 
action constant on frequency is quite weak as compared 
with rapidly decreasing a -1 (u>) as follows from (A3). The 
function ^a^ l {uj) takes values within [—2,0.5] interval 
at the frequencies close to resonant frequency luq. There- 
fore, dispersion equation (10) has a real solution for qa/n 
within a mini-band of frequencies near the resonant fre- 
quency of the scatterers ojq, and this solution is a decay- 
ing function of frequency which corresponds to the back- 
ward wave (the group velocity v g — ^ is negative) . The 
obtained result, of course, confirms existence of backward 
wave below cutoff of the hollow waveguide predicted and 
experimentally demonstrated in [1, 2]. 

In contrast to Fig. 8, the dependencies of normal- 
ized interaction constants C(fc,q) with q = (q,0, 0) T , 
q = (0,7r/a, q) T and q = (q,ir/a,ir/a) T are decaying 
functions of q for the fixed values of k < ir/a. It means 
that solutions of dispersion equations (11), (12) and (13) 
are forward waves for any resonant frequency of the scat- 
terer below ir / (a^/£ Mo)- The Fig. 9 shows dependence of 
the real part of normalized interaction constant C(fe,q) 
with q = (q, 0, 0) T on the normalized frequency ka/ir and 
on the propagation constant qa/ir. The dependencies for 
the cases q = (0, ir/a, q) T and q = (q, Tr/a, ir/a) T are not 
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FIG. 9: Dependence of the real part of normalized interac- 
tion constant C(k, q)a 3 with q = (g, 0,0) T (corresponding to 
longitudinal orientation of magnetic scatterers) on normalized 
frequency ka/n and propagation constant qa/n. 



shown in order to reduce size of the paper. 
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FIG. 10: Dispersion curves for metallic waveguides loaded by 
magnetic scatterers: exact solution (thick line) and effective 
medium filling approximation (dashed line) for transverse ori- 
entation, and exact solution (thin line) for longitudinal orien- 
tation. The effective medium filling model for the longitudinal 
case is not applicable. 

The dispersion curves for the case of magnetic scatter- 
ers are presented in Fig. 10. The thick solid line repre- 
sents the dispersion curve for the transverse mode. It is 
obtained by numerical solution of transcendental disper- 
sion equation (10). The dashed line shows the result pre- 
dicted by the model of effective medium filling (17). The 
significant frequency shift between the exact and approx- 
imate solutions is observed. Also, the effective medium 
model gives wittingly wrong results with q > n/a in the 



region ka < 1.0055, and incorrectly describes group ve- 
locity for q > 7r/(2a) (for example, it does not describe 
the Bragg mode with zero group velocity at the point 
qa = it). The dispersion curve for the longitudinal mode 
obtained by numerical solution of equation (11) is repre- 
sented by thin line in Fig. 10. As it was mentioned above, 
the effective medium filling model can not be applied for 
description of this mode. 




ka 



FIG. 11: Dispersion curves for metallic waveguides loaded 
by electric scatterers: exact solution (thick line) and effective 
medium filling approximation (dashed line) for transverse ori- 
entation, and exact solution (thin line) and effective medium 
filling approximation (thin dashed line) for longitudinal ori- 
entation. 

The dispersion curves for the case of electric scatter- 
ers are presented in Fig. 11. The thick and thin solid 
lines shows dispersion curves for transverse and longi- 
tudinal modes, obtained by numerical solution of dis- 
persion equations (12) and (13), respectively. The dis- 
persion curves provided by effective medium models for 
transverse and longitudinal modes (formulae (20) and 
(21)) are plotted by thick and thin dashed lines, respec- 
tively. The comparison of exact and approximate solu- 
tions shows that the model of effective medium filling 
gives qualitatively right prediction of dispersion curves 
behavior in the case of electric scatterers as well as in the 
case of transverse magnetic scatterers. The drawbacks 
are also the same: wrong group velocity q > -k /(2a) and 
wittingly wrong results with q > n/a at some frequencies. 

The Figs. 10 and 11 demonstrate that the waveguides 
loaded by electric and magnetic resonant scatterers sup- 
port modes within the mini-bands below the cutoff fre- 
quency of the hollow waveguide. The modes are forward 
waves, except the case of transverse magnetic scatterers 
when the mode is forward wave. The bandwidth in the 
case of transverse electric scatterers is of the same order 
with bandwidth for magnetic scatterers with the same 
parameters obtained with help of duality principle, but 
the bandwidths for the longitudinal modes are signifi- 
cantly narrower than for the transverse ones. 
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IV. DISCUSSION 

In the papers [1, 2] the term 'subwavelength waveg- 
uide' was applied to a rectangular waveguide with small 
transversal dimensions as compared to the wavelength 
in free space. However, there is a lot of other works 
[9-17]. in which the term 'subwavelength waveguiding' 
means the propagation of a wave along the chain of elec- 
trically small nearly-resonant particles below the diffrac- 
tion limit. In this case the transversal size of the spatial 
domain, in which the field of the guided mode is concen- 
trated, is much smaller that the wavelength in free space. 
Therefore, this mechanism of the wave transmission is 
considered as prospective for subwavelength imaging. 

Since the field of the mode guided along the chain 
of resonant scatterers is concentrated within a subwave- 
length cross section, the presence of the metal walls even 
at the rather small distance turns out to be not crucial 
for the existence of the guided wave. Thus, any waveg- 
uide periodically loaded by the scatterers can be consid- 
ered as a subwavelength waveguide formed by the chain 
of nearly resonant scatterers, whose dispersive properties 
are perturbed by the metal walls. These walls can be de- 
scribed in terms of the image chains forming an infinite 
lattice. However, the wave propagates along the same 
direction in every image chain, and the transversal wave 
numbers q x = ir/a or q y = ir/b describe the transversal 
phase distribution of the wave propagating along z but 
not the energy transport across z. The main question is 
how the image chains of scatterers interact with the real 
chain, and how this interaction influences its dispersive 
properties. 

For that purpose let us compare Fig. 6 with Figs. 
10 and 11. We can conclude that the presence of metal 
walls around the chain of resonant scatterers produces 
the following effects: 

• It decreases the group velocity and the frequency 
band of the guided mode which corresponds to the 
longitudinal orientation of dipoles. 

• It cancels the two-mode regime for the transverse 
orientation of dipoles, so that the dispersion branch 
becomes backward for magnetic scatterers, and for- 
ward for electric ones. 



Finally, we would like to emphasize, that the width 
of the pass band for the waveguide loaded by transver- 
sal electric scatterers has the same order as in the case 
of transversal magnetic scatterers (Marques waveguide 
[1, 2]) if the loading scatterers have parameters obtained 
using duality principle from each other. So, the loading 
by electric scatterers could be an alternative and even 
more appropriate solution for the waveguide miniatur- 
ization than the design suggested in [2] for this purpose. 
V. CONCLUSION 



The dispersion properties of rectangular waveguides 
loaded by resonant scatterers (magnetic and electric 
ones) have been studied. The waveguide problem has 
been transformed using the image theory into the eigen- 
mode problem of an auxiliary three-dimensional electro- 
magnetic crystal. The dispersion properties of such elec- 
tromagnetic crystal have been modelled using the local 
field approach. It has been revealed that not only mag- 
netic but also electric resonant scatterers allow to ob- 
tain mini pass band below cutoff frequency of the hol- 
low waveguide. The corresponding mini-band turns out 
to be of the same order as for magnetic scatterers with 
same individual frequency dispersion. So, the electric 
scatterers (inductively loaded short wires) could be also 
prospective for the waveguide miniaturization as well as 
split-ring- resonators. It has been shown that the loading 
by scatterers with longitudinal orientation of dipole mo- 
ments also allows to obtain the mini-band of propagation 
below the cutoff frequency of the hollow waveguide, but 
the width of this band is significnatly narrower as com- 
pared to the case of transverse orientation. The observed 
effects arc explained in terms of the subwavelength guid- 
ing properties of the single chains of scatterers. This ex- 
planation is supported by comparison of dispersion prop- 
erties of the loaded waveguides and the chains of the res- 
onant scatterers in free space. Results of our theory are 
in good agreement with the known literature data. For 
the chains of resonant dipoles the results from [12] are 
reproduced. For the rectangular waveguide loaded by 
split-ring-resonators the same result as in [1] has been 
obtained. 
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L is inductance of the ring (we assume that both rings 
have the same inductance): 







NT) 


- 2 







L = [i r 

M is mutual inductance of the two rings: 

'4 s 



M = n r 



(l-01og(| 



2 + C 



w + d 



C r is the effective capacitance of the SRR: 

C r — Eq— arccosh ( ^ ) , 
7T \ d J 

r is the radiation reaction factor: 
„ Aujk 3 



6tt^ 

r is the inner radius of the inner ring, w is the width of 
the rings, d is distance between the edges of the rings (see 
Fig. 2. a), £o and /Jo are permittivity and permeability of 
the host media, and k = uo^e ^ is the wave number 
of the host medium. The presented formulae arc valid 
within the frame of the following approximations: w, d <C 
r and the splits of the rings are large enough compared 
to d. Also, we assume that SRR is formed by ideally 
conducting rings (no dissipation losses). 

The magnetic polarizability (Al) takes into account 
the radiation losses and satisfies to the basic Sipe- 
Kranendonk condition [22, 26, 27] which in the present 
case has the following form: 



APPENDIX A: POLARIZABILITIES OF 
RESONANT SCATTERERS 

1. Split-Ring Resonators 

The SRR considered in [3, 5, 18] is a pair of two copla- 
nar broken metal rings (see Fig.2.b). Since the two loops 
of an SRR are not identical the analytical models of it are 
rather cumbersome [4, 24]. In fact, such SRR can not be 
described as a purely magnetic scatterer, because it ex- 
hibits bianisotropic properties and has resonant electric 
polarizability [4, 24] (see also discussion in [25]). How- 
ever, the electric polarizability and bianisotropy of SRR 
is out of the scope of this paper. We neglect these effects 
and consider an ordinary SRR as a magnetic scatterer. 
The analytical expressions for the magnetic polarizabil- 
ity a{uo) of SRRs with geometry plotted in Fig. 2. a were 
derived and validated in [24]. The final result reads as 
follows: 



Au; 2 



■ j wT ' 



lid 

L + M' 



(Al) 



where u>o is the resonant frequency of magnetic polariz- 
ability: 

1 



UJn = 



(L + M)C r : 



Imja- 1 ^)} 



6-7T/XO 



(A2) 



In our analysis we operate with the inverse polarizabil- 
ity a~ 1 (oj), thus, we rewrite (Al) in the following form: 



1 +j 



67T^ 



(A3) 



2. Inductively Loaded Short Wires 

An inductively loaded short wire is shown in Fig. 2.b. 
The electric polarizability a e of an inductively loaded 
wire following the known model [19] has the form: 



3 



l-u 2 /u 2 
A-u 2 /u 2 



6tTEq 



(A4) 



y wirc = kIeq/ log(2l/ro) is the capacitance of the 
wire, uj Q = ^LCv, 



where C, Viii 

is the resonant frequency, L is the 
inductance of the load, I is the half length of the wire 
and ro is the wire radius. 

It is clear, that at the frequencies near the resonance 
the polarizability of LSW has the form 



a-V) 



UJ 2 



+ 3 



67T£o ' 



(A5) 
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with A e = ? 2 C w i ro , which is similar to (A3). Moreover, if 
Ae/£o = A/(1q then using duality principle the magnetic 
dipole with polarizability a (A3) can be transformed to 
the electric dipole with polarizability a e (A3), and vice 
versa. This means that it is enough to consider only one 
type of resonant scatterers. In the present paper we have 
chosen magnetic ones to be principal. The case of electric 
scatterers was obtained using the duality principle with 

A = fX A e /e . 

APPENDIX B: INTERACTION CONSTANTS OF 
THE CHAINS 



where 

1+ = l _ e -](k+q)a t~ — I — e -j{k-q)a 

t + = \- e j(k+q)a , t- = 1 - e J (k - q)a . 

As to C y , the series of the order 1/m in (B2) can be 
obtained in the closed form using the tabulated formula 
(see [21], Appendix): 

_ = -lDg(l- e -") (B4) 



The initial expressions for the interaction constants 
C x , y entering (7) follow from (4) and (5) and read as 
follows: 



E l + jka\ 
27ra 3 |m 



m \ e -j(k\m\+qm)a 



(Bl) 



+ oo 



= — y 

7ra 3 ^ 



m—l 



in 



jka 



e -jkam cos ( gam); 



for the longitudinal polarization, and 
k 2 a 2 m 2 — jka\m\ — 1 



= 2W ^ 

m—l 



47ra 3 |TO| 3 



j(k\m\+qm)a 



(B2) 



k 2 a 2 



m 



jka 
m 2 



m° 



-jkam 



cos(qam) 



for the transverse one (see e.g. [12]). 

Note, that C x includes only near-field terms (of the 
order 1/R 2 and 1/-R 3 ). In contrast to C x , the transverse 
interaction constant C y includes also the wave terms (of 
the order 1/R) which corresponds to the slowly converg- 
ing series. It makes the direct numerical summation of 
(B2) to be not efficient. The series in (Bl) have better 
convergence, but it is also not enough for rapid calcula- 
tions. 

The application of acceleration technique done in [20] 
offers a more cumbersome expression for C x {k, q, a) than 
(Bl), but it is better for numerical calculations since the 
series converges very rapidly: 



47TO 3 



(2jka + 3)m + 2 
' m 3 (m + l)(m + 2)' 



-jkam cos ( qam ) 



-(jka + 1) (4 logt+ + t 2 _ log* - + 2e lka cos(qa)) (B3) 



-2jka (t+ \ogt + + t- logi") + {7jka + 3) 



log 



2sin^ 



7T — 7 



where 7' = 27r{7/(27r)} and we use notation {x} for frac- 
tional part of variable x. The rest part of the expression 
(B2) is simply proportional to C x . Thus, C y can be eval- 
uated as follows: 



C y (k, q, a) — —- — log 1 2 (cosfca — cosqa)| 



k 2 , 



(k + q)a 



(k — q)a 



2tt J { 2ir 
In the works [17, 20] it was shown that 



C x (k,q,a)/2. 

(B5) 

k 2 ) , (B6) 



where 



\q m \<k 



27rm 



q - 



From the formula (B2) using some algebra it follows 
that 



(B7) 



\q m \<k 



The expressions (B6) and (B7) demonstrate energy 
transformations happening in the chains. A single scat- 
terer radiates cylindrical wave and that is why its polar- 
izability has radiation losses which can be described by 
Sipe-Kronendonk condition (A2). Being arranged into 
the regular chains the scatterers acquire effective polar- 
izability (with respect to the external field) of the form: 



a x,y — \_Oi 1 H C x . y ~\ 



(B8) 



where indices x and y correspond to longitudinal and 
transverse orientations of scatterers in the chain, respec- 
tively. In accordance to (B6) and (B7) one can formulate 
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the following analogues of Sipe-Kroncndonk condition for 
effective polarizabilitics of the scatterers in the chains: 



Imjc^ 1 } 



1 



E ( fc2 -^) 



\q m \<k 



(B9) 



formula deducted in [20] : 



C(k, q, a, 6, c) = - ^ —K a (p m bri) cos(q y bn) 

n=l Rc(p m )^0 



Im{a/} = 



8/ioa 



E fe + fc2 )- 

!<? m |<fc 



(BIO) 



+00 +00 9 _o'fc( mn ),- 

e 3K * - cosg z c 



+ E E 



n ^_ GO 2jabkz cos kz 'c— cos q z c 



(Cl) 



Note, that terms fc 3 / (67r) are canceled. It is clear, that 
if there are no such index m that \q m \ < k (like it happens 
for example if k < q < lix j a — k for k < ir/a) then ef- 
fective polarizabilities turn out to be purely real and the 
chain itself does not radiate. This regime corresponds to 
the case of guiding modes and it makes dispersion equa- 
tion (7) real valued one. If there are some indices m that 
\q m \ < k then the effective polarizabilites acquire non- 
zero imaginary part which give evidence that the chain 
radiates cylindrical waves. The number of such waves 
corresponds to the number of indices m fulfilling to the 
relation \q m \ < k. If \q\ < k < ir/a then the chain ra- 
diates only one cylindrical wave which can be treated as 
the main diffraction lobe of this periodical array. For 
the higher frequencies the grating lobes appears and all 
of them make its contribution to expressions (B9) and 
(BIO). 



APPENDIX C: INTERACTION CONSTANT OF 
AN ORTHORHOMBIC LATTICE 

For effective numerical calculation of interaction con- 
stant C(k, q, a, b, c) defined by (9) we use the following 



r, 2 I 1 + °° 
V- Prn 1 + y 

Re( Pm )=0 \J Kz n=l 



b l m b 



.3 1 



7T7J 87T 3 n 3 



+C x (k,q x ,a), 
where C x is given by (B3) and 

Urn) _ , u(n) 



2Tm 



(k^y 



fc 2 , 



b ' 



= -J\l(k { X m) ) 2 +(ki n) ) 2 -k 2 . 



